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The evaluation of the hadronic contribution to the muon magnetic anomaly a M is revisited, taking advantage of 
new experimental data on e + e~ annihilation into hadrons: SND and CMD-2 for the tt + tt~ channel, and BABAR 
for multihadron final states. Discrepancies are observed between KLOE and CMD-2/SND data, preventing one 
from averaging all the e + e~ results. The long-standing disagreement between spectral functions obtained from r 
decays and e + e~ annihilation is still present, and not accounted by isospin- breaking corrections, for which new 
estimates have been presented. The updated Standard Model value for a M based on e + e~ annihilation data is 
now reaching a precision better than experiment, and it disagrees with the direct measurement from BNL at the 
3.3a level, while the T-based estimate is in much better agreement. The r/e + e~ discrepancy, best revealed when 
comparing the measured branching fraction for r~ — ■» ty~tt°u t to its prediction from the isospin-breaking-corrected 
e + e~ spectral function, remains a serious problem to be understood. 



1. Introduction 

Hadronic vacuum polarization (HVP) in the 
photon propagator plays an important role in 
many precision tests of the Standard Model. This 
is the case for the muon magnetic anomaly = 
(g fl — 2)/2, where the HVP component is the lead- 
ing contributor to the uncertainty of the Stan- 
dard Model prediction. The HVP contribution is 
computed by means of a dispersion relation as an 
integral over experimentally determined spectral 
functions. It is a property of this dispersion re- 
lation that the ttit spectral function provides the 
major part of the total HVP contribution, so that 
the experimental effort focuses on this channel. 

Spectral functions are directly obtained from 
the measured cross sections of e + e _ annihilation 
into hadrons. The accuracy of the HVP predic- 
tions has therefore followed the progress in the 
quality of the data [ Q] it relies on. Because 
the data quality was not always suitable, it was 
deemed necessary to resort to other sources of in- 
formation. One such possibility was the use of 
the vector spectral functions [ 2J derived from 
the study of hadronic r decays [ [3] for the en- 
ergy range less than m T ~ 1.8GeV/c 2 . For this 
purpose, the isospin rotation that leads from the 



charged r to the neutral e + e~ final state has to be 
thoroughly corrected for isospin-breaking effects. 

Also, it was demonstrated that essentially per- 
turbative QCD could be applied to energy scales 
as low as 1-2 GeV [ IU [5] , thus offering a way to 
replace poor e + e _ data in some energy regions 
by a reliable and precise theoretical prescription [ 

El El El EH EH- 

Detailed reanalyses including all available ex- 
perimental data have been published in Refs. [ 
[T2l [T3l ITU UM ITS] , taking advantage of precise 
results in the 7T7t channel from the CMD-2 ex- 
periment [ [T7| and from the ALEPH analysis of 
r decays [ [TB] , and benefiting from a more com- 
plete treatment of isospin-breaking corrections [ 
HH1I5D]. With the increased accuracy of the e + e~ 
data a discrepancy with the isospin-breaking- 
corrected t spectral functions was found [ IT2"] . 
thus leading to inconsistent predictions for the 
lowest-order hadronic contribution to a^. The 
dominant contribution to this discrepancy stems 
from the 7r + 7r~ channel, although another dis- 
crepancy occurs in the tt + tt~2tt° mode. 

Improvements in the HVP calculation are 
needed in order to match the present experimen- 
tal accuracy on a M from the BNL experiment [ 
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£5, 

a M = (11 659 208.0 ±6.3) x 10~ 10 . (1) 

In this paper I revisit the input to the HVP dis- 
persion integral in the light of new experimental 
data on e+e" -> tt+tt" from SND [[22] and CMD- 
2 [ E3] EH [25] , and on multihadron final states 
from BABAR [ [261 EH EH [29] using the radiative 
return technique [[3D]- These new measurements 
represent a significant step forward, as they over- 
come in precision previous determinations in the 
same channels. 

2. Muon Magnetic Anomaly 

It is convenient to separate the Standard Model 
(SM) prediction for the anomalous magnetic mo- 
ment of the muon into its different contributions, 

^ M = a^+a^ + ajr k , (2) 
with 

„had had.LO , had, HO , had.LBL /o\ 

% = a M + % + % ' , (3) 

and where a® ED = (11 658 471.9 ± 0.1) x 10" 10 
is the pure electromagnetic contribution [ 31j, 
a had,LO j g foe lowest-order HVP contribution, 
a had,HO = (_9.8 ±0.1) x 10" 10 is the correspond- 
ing higher-order part [ El El [14], and a™ cak = 
(15.4 ± 0.1 ± 0.2) x 10" 10 , where the first er- 
ror is the hadronic uncertainty and the second is 
due to the Higgs mass range, accounts for correc- 
tions due to exchange of the weakly interacting 
bosons up to two loops [ [33] . For the light-by- 
light (LBL) scattering part, a^ ad ' LBL , we use the 
value (12.0 ± 3.5) x 10~ 10 from the latest evalu- 
ation [ 34, 35 , slightly corrected for the missing 
contribution from (mainly) the pion box. 

Owing to unitarity and to the analyticity of 
the vacuum-polarization function, the lowest or- 
der HVP contribution to can be computed via 
the dispersion integral [ 136] 

oc 

«r LO = ^ jds^i^Hs), ( 4) 

where K(s) is a well-known QED kernel, and 
R>°'(s) denotes the ratio of the "bare" cross 



section for e + e~ annihilation into hadrons to 
the pointlike muon-pair cross section. The 
bare cross section is defined as the measured 
cross section corrected for initial-state radiation, 
electron-vertex loop contributions and vacuum- 
polarization effects in the photon propagator. 
However, photon radiation in the final state is in- 
cluded in the bare cross section defined here. The 
reason for using the bare (i.e., lowest order) cross 
section is that a full treatment of higher orders 
is anyhow needed at the level of a M , so that the 
use of the "dressed" cross section would entail the 
risk of double-counting some of the higher-order 
contributions. 

The function K(s) ~ 1/s in Eq. Q gives a 
strong weight to the low-energy part of the in- 
tegral. About 91% of the total contribution to 
a had,LO j s accumulated at center-of-mass energies 
y/s below 1.8 GeV and 73% of a ^ LO is covered 
by the irn final state, which is dominated by the 
p(770) resonance. 

3. The 2vr Input Data 

A detailed compilation of all the experimen- 
tal data used in the evaluation of the disper- 
sion integral ([4]) prior to 2004 is provided in 
Refs. [ [12] . Also discussed therein is the 
corrective treatment of radiative effects applied 
to some of the measurements. The r 2tt spec- 
tral function is obtained by averaging the re- 
sults from ALEPH [0], CLEO [[37] and OPAL [ 
138] . which exhibit satisfactory mutual agreement. 
Since 2004, new cross section measurements be- 
came available from KLOE [ [39] using radiative 
return at DAPHNE and from the annihilation ex- 
periments SND [ E2] and CMD-2 [ E3 E2 E5] at 
Novosibirsk. 

The comparison between all new e + e~ results 
and the combined r spectral function is given in 

Fig. m 

A few remarks are in order: 

• Revision of the radiative corrections applied 
the CMD-2 (94-95 data) and SND data led 
to corrections amounting up to 3%. 

• The revised SND and CMD-2 (also includ- 
ing new data released in 2006) spectral 
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Figure 1. Relative comparison of the tt + tt spectral functions from e + e -annihilation data and isospin- 
breaking- corrected r data, expressed as a ratio to the t spectral function. The shaded band indicates 
the errors of the r data. The e+e~ data are from KLOE [W$, CMD-2 [\T^, CMD, OLYA and DM1 
(references given in Ref. [\12lj). The right hand plot emphasizes the region of the p peak. 



functions agree within errors. It should 
be pointed out that both analyses now use 
the same radiative correction package, in- 
troducing a full correlation between the two 
data sets. 

• The high-statistics KLOE data do not agree 
with SND and CMD-2, mostly through a 
discrepancy in the p lincshape: KLOE is 
higher below the peak and becomes lower 
above. 

• A significant discrepancy, most visible 
above the p peak, but present almost ev- 
erywhere, is found between r and the e + e~ 
data. 

Considering different p~ and p° masses as an 
additional isospin-breaking correction of the r 
spectral function would improve the comparison 
in the p resonance peak, but in that case the 
discrepancy should fade away for masses above 
m p + r p /2, which is not observed. 

Finally, some improvement of the isospin- 
breaking corrections has been proposed in Ref. [ 
140] , A contribution from the pcun vertex with 
lo — ► 7r°7 was not included in the treatment used 



so far [ [20] , as it occurs from higher order in Chi- 
ral Perturbation Theory. It was however found to 
be significant, but its effect amounts to only 20% 
of the observed discrepancy. 

4. Testing CVC 

Measurement of branching fractions in r de- 
cays are more robust than the spectral functions, 
as the latter ones depend on the experimental res- 
olution and require a numerically delicate unfold- 
ing. It is possible to relate the measured branch- 
ing ratios for r~ — > V~v T , where V is any vector 
final state, to their expectations from CVC us- 
ing e + e~ spectral functions, duly corrected for 
isospin breaking. In this way we do not anymore 
rely on the shape of the r spectral function and 
instead concentrate on the relative normalization 
and the isospin-breaking corrections. 

The result of the test for the tt~it channel is 
shown in Fig. O It shows a large discrepancy be- 
tween the average r branching ratio and the CVC 
prediction. The difference [B T — Bcvc]^" = 
(0.92 ± 0.21)% is 4.5cr away from zero. In rel- 
ative terms, the discrepancy is a 3.6% effect, 
about twice the already applied isospin-breaking 
correction, dominated by (expected to be) well- 
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Figure 2. TTie measured branching ratios for 
t~ — > v T ix~ifi compared to the prediction from 
the e + eT — > tt + it~ spectral function applying the 
is o spin-breaking correction factors discussed in 
Ref. [ \13f. The measured branching ratios are 
from ALEPH [U3j, CLEO [ JTj and OPAL [ 
The L3 and OPAL results are obtained from 
their hir° branching ratio, reduced by the small 
Kir a contribution measured by ALEPH [ 4-31 and 
CLEO [\44$. 



controlled short-distance effects. 

The same exercise can be applied to other chan- 
nels. For the 7r37r° final state, the difference is 
consistent with zero, — (0.08±0.11)%. The result 
for the 3mr° mode is bad, (0.91±0.25)%, amount- 
ing to a relative 20% discrepancy, certainly much 
beyond any reasonable isospin-breaking effect. 
Here the quality of the data on e + e~ — > tt + tt~2tt° 
is in doubt, as they show a considerable spread 
between the experiments. 

5. New Multihadron Data 

Results from the BaBAR [ HH HO ex- 
periment are being produced on the multihadron 
final states using radiative return [ [30] . They are 
part of a program designed to cover most exclu- 
sive annihilation processes in the few GeV energy 
range, taking advantage of the large initial ccntre- 
of-mass energy of 10.6 GeV. Hard-radiated pho- 
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Figure 3. The measured cross section for e + e — ► 
7r + 7T~7r from BABAR [ (indicated by full 
sguares) compared to previous measurements (see 
references in Ref. [\13f). 



tons are detected at large angle, together with 
the hadronic system byproducts, so that the full 
final state can be kinematically constrained (see 
the discussion in Ref. [ 29]). These results signif- 
icantly improve the corresponding contributions 
to a^ ad , since earlier data in the 1.4-3 GeV en- 
ergy range were of poor statistical and systematic 
quality. 

The BABAR results are compared to other data 
in Figs. for the e+e~ -> tt+tt-tt , 2tt + 2tt-, 
37r + 37r - , and 2tt + 2tt^2tt° cross sections. Besides 
a generally good agreement with previous exper- 
iments, some important differences are seen for 
7T + 7r~7r° and the 6-pion states with the results 
obtained at DCI. The impact of these new mea- 
surements is quantified in Table [1] 

6. Results 

During the previous evaluations of a^ ad ' LO , the 
results using respectively the r and e + e~ data 
were quoted individually, but on the same footing 
since the e + e~-based evaluation was dominated 
by the data from a single experiment (CMD- 
2). The confirmation of the r/e + e~ discrepancy 
by SND and KLOE may suggest to prefer the 
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Figure 4. The measured cross section for e + e — > 
27T+27T - /rom BaBAR /[f^ (indicated by full cir- 
cles) compared to previous measurements and re- 
sults from t~ — > 7T~37r° (see references in Ref. [ 
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Figure 6. The measured cross section for e 
27T + 27r-27r° from BABAR [ 28} (indicated by tri- 
angles pointing below) compared to previous mea- 
surements (see references in Ref. [Mttf). 
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Figure 5. The measured cross section for e + e~ — » 
37T + 37r _ from BABAR [ 281 (indicated by open cir- 
cles) compared to previous measurements (see ref- 
erences in Ref. [HE/)- 



Table 1 

The contribution of some multipion processes to 
a^ ad integrated from threshold to 1.8 GeV for the 
older experiments (references in Ref. and 
including the new BABAR ISR results \2~T\ 
\28\ \ 29$ . Values are given in units o/10 -10 . 



process 



older exp. with BABAR 



7T+7T 7T° 

2ir + 2n- 

37T+37T - 

27T+27r^27r 



2.45 ±0.26 
14.20 ±0.90 
0.10 ±0.10 
1.42 ±0.30 



3.25 ±0.09 
13.09 ±0.44 
0.11 ±0.02 
0.89 ±0.09 
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e + e~ -based result until a better understanding 
of the dynamical origin of the observed effect is 
achieved. This discrepancy is a challenging prob- 
lem, which may itself turn out to be of fundamen- 
tal importance. The present update is therefore 
only based on e + e~ data, dominated by the latest 
results from CMD-2 and SND which are in good 
agreement. The KLOE data show a systematic 
trend which is not explained at the moment, so we 
do not use them as the resulting increase in pre- 
cision would not be trustworthy. Further studies 
by KLOE are ongoing, in particular a determina- 
tion of the 7T7r//i^ ratio where several systematic 
effects cancel, which should lead a significant im- 
provement of the systematic uncertainty [ [JS] . 

The preliminary estimate of the integral (|4|) 
given below includes one additional improvement 
with respect to Ref. [02]: perturbative QCD is 
used instead of experimental data in the region 
between 1.8 and 3.7 GeV, where non-perturbative 
contributions to integrals over differently weighed 
spectral functions were found to be small [ [7J. 
This results in a reduction of a^ ad,LO by — 1 x 
10 -10 . All contributions to the dispersion inte- 
gral where no new input data are available are 
taken from Ref. [ [T^] . 

The R values from data and QCD are displayed 
in Fig. but not yet updated with the BABAR 
multipion data. For masses larger than 1.8 GeV, 
except in the cc threshold region from 3.7 to 5 
GeV, the QCD prediction is used. Agreement be- 
tween QCD and data is good. The contributions 
of the different exclusive channels and of the con- 
tinuum are given in Table [2] 

The e + e~-based result for the lowest order 
hadronic contribution is 

a had,LO = ( 6 90.8±3.9±1.9 rad ±0.7Q CD )xl0- 10 , (5) 

where the second error is due to our treatment 
of (potentially) missing radiative corrections in 
the older data [[13]. For comparison, the r-based 
result [ 1 1 2 j can be updated using the new e + e~ 
BABAR data for the channels other than 2tt or 4ir 
yielding the value (710.3 ± 5.2) x 1CT 10 . Adding 
to the e + e~ result (0 the QED, higher-order 
hadronic, light-by-light scattering, and weak con- 
tributions given in Section [51 one finds 

al M = (11 659 180.5 ±4.4 had , LO+ HO 
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Figure 8. Comparison of the result {6]) with the 
BNL measurement [\21}j. Also given are our pre- 
vious estimates where the triangle with the 
dotted error bar indicates the T-based result, as 
well as the estimates from Refs. [ \14\ \16f , not 
yet including the KLOE data. 

± 3.5had,LBL ±0.2 QED+EW ) x 10- 10 . (6) 

This value can be compared to the present mea- 
surement ([1]); adding all errors in quadrature, the 
difference between experiment and theory is 

a™ P -a™ = (27.5 ±8.4) x 10~ 10 , (7) 

which corresponds to 3.3 "standard deviations" 
(to be interpreted with care due to the dominance 
of systematic errors in the SM prediction). A 
graphical comparison of the result ^ with pre- 
vious evaluations and the experimental value is 
given in Fig. [8] 

7. Conclusion and Perspectives 

In spite of the new and precise data on the two- 
pion spectral function from CMD-2 and SND, and 
on multihadron cross sections from BaBAR, the 
lowest order hadronic vacuum-polarization con- 
tribution remains the most critical component in 
the Standard Model prediction of a^. Yet, for 

1 Results similar to ours have just appeared in Ref. [I46|. 



7 




• Crystal Ball o MARK J 



5 6 7 8 9 10 11 12 13 14 

Vs (GeV) 



Figure 7. The R values from data and QCD used in the dispersion integral. The shaded band below 
2.1 GeV represents the sum of the exclusive channels, not yet including the new BABAR results. Narrow 
resonances are indicated as dashed lines. All data points shown correspond to inclusive measurements. 
The cross-hatched band gives the prediction from QCD used in the calculation from 1.8 to 3.7 GeV and 
above 5 GeV. 



the first time in recent years, the accuracy of the 
prediction exceeds that of the experiment. One 
should not forget however that the theoretical er- 
ror is completely propagated from the systematic 
uncertainties of the input experiments, the es- 
timation of which we totally depend. Also the 
evaluation of the systematic uncertainty on the 
hadronic light-by-light contribution is more sub- 
ject to caution. 

The discrepancy between the 2tt spectral func- 
tions obtained from r decays and from e + e~ an- 
nihilation is still unresolved: it affects both the 
overall normalization and the shape. A partic- 
ularly important test of the relative normaliza- 
tion is obtained, when comparing the measured 
branching fraction for r~ — > 7r _ 7r°i/ r to its Stan- 



dard Model prediction from the e + e~ spectral 
function, corrected for small isospin-breaking ef- 
fects. Here the discrepancy amounts to 4.5 stan- 
dard deviations. 

In view of this problem and the fact that both 
Novosibirsk experiments agree at the 1% level, 
the hadronic contribution is computed only with 
e + e~ data, excluding for the moment KLOE un- 
til their disagreement with Novosibirsk is under- 
stood. The choice of using only two experimental 
inputs (CMD-2 and SND) among the four avail- 
able (+KLOE and r decays) is clearly not satis- 
factory and even debatable, but represents in our 
view the most reasonable solution, all arguments 
considered. In this scenario we find that the Stan- 
dard Model prediction of differs from the ex- 
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Table 2 

The contributions of different final states in specified energy ranges to a„ , given in units o/10" 10 . Major 
improvements since Ref. [ \12]j are obtained for (1) the tt + tt~ channel from CMD-2 (here preliminary 
results are used; final data have been published [ 23, 21, \25\l since) and SND f\22f, while KLOE results are 
not included (see text), and (2) the 2tt + 2tt~ and other exclusive channels from BaBAR [Wb\ \W7\ [2~3\. I^Pf . 
The uncertainty for missing radiative corrections, labeled 'rad', only concerns our ad hoc treatment 
applied to older experimental data. 



Modes 


Energy range (GeV) 


~had 
U M 




2771^-0.5 


55.6 ± 0.8 ±0.1 rad 


7T + 7T~ 


0.5-1.8 


449.0 ± 3.0 ± 0.9 rad 


27T+27T- 


2m,r-1.8 


13.1 ± 0.4 ±0.0 rad 




2777^-1.8 


16.8 ± 1.3 ±0.2 rad 


w 


0.3-0.81 


38.0 ± 1.0 ±0.3 rad 




1.0-1.055 


35.7 ± 0.8 ±0.2 rad 


other exclusive 


2777^-1.8 


24.3 ± 1.3 ±0.2 rad 


JMV(2S) 




7.4±0.4±0.0 rad 


R(QCD) 


1.8-3.7 


33.9±0.5 Q cd 


R(data) 


3.7-5.0 


7.2±0.3±0.0 rad 


R(QCD) 


5.0-oc. 


9.9±0.2 QCD 


sum 


2777^-00 


690.8 ± 3.9 ± 1.9 rad ± 0.7 QCD 



perimental value by 3.3 standard deviations. 

We are looking forward to the forthcoming 
results on the two-pion spectral function from 
KLOE and BABAR using the tttt/ pLfx ratio. Since in 
this way vacuum polarization cancels, these data 
will help to reduce the systematic uncertainty due 
to the corrective treatment of radiative effects, 
always problematic in previous experiments nor- 
malized by luminosity. More data from BABAR on 
multihadron final states are also expected soon. 
With new experimental input to the vacuum po- 
larization integrals, the quality of the prediction 
will improve, opening the way to a more precise 
direct determination of [ 27] . Unfortunately, 
some recent prospective work in the US [|48] does 
not cast a bright future in this direction. 
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